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An  elementary  model  of  how  a turbulent  boundary  layer  (b.l.)  generates  noise 
is  analyzed.  Two  similarity  hypotheses  are  given  which  relate  b.l.  turbulence 
to  flow  noise. 

Dynamic  pressures  come  from  the  forces  created  by  the  interaction  of 
eddies  on  the  viscous  sublayer  where  velocity  components  normal  to  the 
wall  are  brought  to  rest.  6>^oi. 

The  forces  are  proportional  to  the  second  power  of  the  eddy  velocities 
and  the  distribution  of  separation  distance  between  pressure  points  on 


EDITION  OF  I NOV  «3  IS  OBSOLETE 


78  07  12  0^Q& 


UNCLASSIFIED 

ECURITV  CLASSIFICATION  OF  THIS  PAGE  Cllhwi  Dmtm  Eni.m* 


SECUniTY  CLASSIFICATION  OF  THIS  PAGEfMTi^n  DmtM  Enfered) 


20.  Abstract  (Continued) 


the  wall  is  the  same  as  the  distribution  of  eddy  sizes. 

It  is  concluded  that  flow  noise,  wall  pressure  f luctutations  and  b.l. 
turbulence  all  have  the  same  spectral  shape  in  frequency  space.  Data  on 
wall  pressures  fluctuations  generally  have  a(Y~ll7^form  due  to  spatial 
averaging  by  the  pressure  sensor.  Corrections  for  pressure  sensor 
averaging  are  given  and  experimental  data  is  given  to  demonstrate  the 
validity  of  the  correction.  / 
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An  elementary  model  of  how  a turbulent  boundary 
layer  (b.l.)  generates  noise  is  analyzed.  Two 
similarity  hypotheses  are  given  which  relate  b.l. 
turbulence  to  flow  noise. 

I.  Dynamic  pressures  come  from  the 

forces  created  by  tae  interaction  of  eddies 
on  the  viscous  sublayer  where  velocity  com- 
ponents normal  to  the  wall  are  brought  to 
rest. 

II.  The  forces  are  proportional  to  the  second  power  of 
the  eddy  velocities  and  the  distribution  of  sepa- 
ration distance  between  pressure  points  on 
the  wall  is  the  same  as  the  distribution  of 
eddy  sizes. 

It  is  concluded  that  flow  noise,  wall  pressure  fluc- 
tuations and  b.l.  turbulence  all  have  the  same 
spectral  shape  in  frequency  space.  Data  on  wall 
pressure  fluctuations  generally  have  a form 

due  to  spatial  averaging  by  the  pressure  sensor. 
Corrections  for  pressure  sensor  averaging  are  given 
and  experimental  data  is  given  to  demonstrate  the 
validity  of  the  correction. 
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NOMENCLATURE 

A area 

c speed  of  sound 

D diameter 

£ frequency,  Hz 

k wave  number,  k = U/2Trf 

L length 

p pressure 

r radius/ sensor  radius 

u,v  velocity  deviation  from  mean  (one  component) 

U mean  velocity 

a coefficient  for  turbulence  inertial  subrange  (~0.55,  1-D) 

6 boundary  layer  thickness 

e rate  of  disipation  of  kinetic  energy  per  unit  volume  of  fluid 

n Kolmogoroff  length,  r,  = (v^/£)^^^ 

X wave  length  X*U/f 

)j  10  ^ metres 

V kinematic  viscosity 

p density 

^ one  dimensional  velocity  energy  spectrum 

b>  radians/sensor  radius 

Subscripts 
a acoustic 

K dlmenslonalized  as  suggested  by  Kolmogoroff 

m measured 

V velocity 
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Introduction 

Flow  noise  is  generated  by  a turbulent  boundary  layer's  (b,l.)  Inter- 
action with  the  rigid  wall.  experimenal  data  reported  by  Haddle  and 

Skudrzyk  (1969)  show  the  flow  noise  falling  off  at  the  rate  of  5 dB/octave. 

The  data  of  Clay  (1973)  illustrates  that  the  spectrum  of  turbulence  falls 
off  at  a rate  of  5 dB/octave. 

The  spectral  properties  of  wall  pressure  fluctuations  are  not  clearly 
understood.  Foxwell  and  Grandshaw  (1968)  report  a correlation  of  f ^ with  an 
f at  the  highest  frequencies  and  Foxwell  (1966)  reports  an  additional  f~^ 
correlation  at  the  lowest  frequencies.  Haddle  and  Skudrzyk  (1969)  indicate 
that  the  flow  noise  is  dominated  by  turbulence  at  frequencies  above  200-500 
Hz. 

A review  of  the  literature  on  wall  pressure  fluctuations  by  the  author 
revealed  much  experimental  data  of  high  quality  but  a lack  of  understanding 
of  the  basic  phenomena  involved.  It  is  the  purpose  of  this  paper  to  discuss 
how  turbulence  generates  flow  noise  and  the  statistically  averaged  spectral 
properties  of  turbulence,  wall  pressure  and  flow  noise.  A correction  for 
spatial  averaging  by  the  sensor  will  be  presented  along  with  a new  inter- 
pretation of  the  measured  wall  pressure  spectra. 

Theory 

A turbulent  b.l.  is  composed  of  a rigid  wall  and  a layer  of  turbulent 
fluid  separated  by  a very  thin  layer  of  fluid  called  the  viscous  sub  layer. 

The  turbulence  is  characterized  by  a Reynold's  stress,  or,  three  dimensional 
(3-D)  vorticies.  One  can  visualize  a single  eddy  as  being  a cup  of  coffee 
stirred  about  the  cylindrical  surface.  To  construct  a model  of  a b.l.  one 
metre  thick,  suppose  one  had  cylinders  of  random  sizes  ranging  from  one  metre 
to  1 mm  in  diameter.  If  one  stacked  them  on  the  floor  on  a piece  of  carpet 
to  make  a wall  about  1 m high  one  has  a simple  representation  of  a cross 
section  of  a turbulent  b.l.  The  floor  represents  the  rigid  boundary,  the 
carpet  the  viscous  sub  layer  and  the  cylinder  the  vortical  motion  of  the 
eddies.  The  actual  eddies  are  3-D  and  therefore  have  simultaneous  vortical 
motion  about  three  orthogonal  axis.  For  this  reason,  actual  eddies  tend  to  be 
spherical  in  shape.  In  the  hypothetical  model,  an  eddy  represented  by  a "drink- 
ing straw"  may  be  a "mixing  bowl"  sized  eddy  in  a cross  section  displaced  a 


UNCLASSIFIED 


UNCLASSIFIED 


- 6 - 


March  16,  1978 
JPC ;pam 


few  cm  from  the  model  cross  section.  3-D  eddies  are  difficult  to  visualize 
so  we  will  use  the  cylinders.  If  one  now  visualizes  a fork  lift  rapidly  push- 
ing on  one  end  of  the  "wall"  of  cylinders  such  that  the  top  cylinders  are 
moving  faster  than  the  bottom  ones;  one  has  completed  the  modeling  of  the 
turbulent  b.l. 

Pressure  fluctuations  are  generated  by  the  interaction  of  the  eddies  and 

the  viscous  sub  layer.  If  the  local  velocity  of  an  eddy,  normal  to  the  wall, 

is  quickly  brought  to  rest  at  the  wall,  a local  force  equal  to  the  total  head 

2 

pressure  of  the  eddy- (1/2  pv  ) is  made  on  the  viscous  sub  layer.  The  energy 
available  to  be  converted  to  acoustic  noise  is  the  kinetic  energy  of  the  eddy. 
The  eddies  on  the  viscous  sub  layer  are  not  circular,  but  are  flattened  where 
they  interact  with  the  viscous  sub  layer.  Where  the  eddies  are  flattened 
there  is  a local  high  and  low  pressure.  The  local  "high"  pressure  areas  are 
randomly  spaced  on  the  wall  based  on  the  random  size  and  packing  of  the  eddies. 
The  amplitude  of  the  local  pressures  are  random  also.  The  turbulent  flow  is  a 
stochastic  process  but  it  has  clearly  defined  statistical  properties.  There 
are  Interactions  between  eddies  also  but  the  main  interaction  is  between  eddies 
and  the  viscous  sub  layer. 

A local  "high  pressure"  area  immediately  acts  to  relax  the  system  to  make 
the  pressure  uniform.  The  relaxation  is  performed  by  moving  fluid  away  from 
the  high  pressure  area  through  force  equals  mass  times  acceleration.  The 
acceleration  process  travels  "spherically"  outward  at  the  speed  of  sound  and 
is  called  a sound  wave.  The  collection  of  a large  number  of  these  random 
sound  waves  is  called  "flow  noise." 

It  is  proposed  that  two  similarity  hypotheses  relate  flow  noise  to 
turbulence: 

I.  Dynamic  pressures  come  from  the  forces  created  by  the 

Interaction  of  eddies  on  the  viscous  sub  layer  where  velocity 
components  normal  to  the  wall  are  brought  to  rest. 

II.  The  impulsive  forces  are  proportional  to  the  square  of  the 
eddy  vortical  velocities  and  the  distribution  of  separation 
distance  between  pressure  points  on  the  wall  is  the  same  as 
the  distribution  of  the  eddy  sizes. 
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One  consequence  of  these  hypotheses  is  that  the  power  spectra  of  the 

b.  1.  velocity,  wall  pressure  and  flow  noise  plotted  as  magnitude  vs 

frequency  will  have  a single  spectral  form.  In  this  paper,  flow  noise  is 

defined  as  acoustic  noise  originating  from  a turbulent  b.l.  These  acoustic 

signals  are  superimposed  on  the  flow  noise  to  make  the  total  radiated  noise 

field.  The  wave  lengths  of  turbulence  and  flow  noise  are  not  the  same.  The 

wave  length  of  turbulence  and  wall  pressure  is  = U/f  while  the  wave  length 

of  flow  noise  is  X = c/f.  The  theory  of  turbulence  discussed  below  assumes 
a 

incompressible  flow  so  U/c  «1. 

Spectral  Form 

The  spectral  form  of  turbulence  is  well  established  both  theoretically 
and  experimentally.  Kolmogoroff  (1941)  proposed  two  hypotheses  on  th*' 
similarity  of  all  turbulent  velocity  fields.  The  hypotheses  assume  that  the 
flow  fields  are  incompressible,  isotropic,  homogeneous  and  Re  -►<».  The  first 
hypothesis  states  that  the  distribution  of  velocity  fluctuations  is  uniquely 
determined  by  the  kinematic  viscosity  (v)  and  the  mean  rate  of  energy  dissi- 
pation per  unit  mass  (e).  This  hypothesis  implies  that  fluid  motion  is 
dominated  by  viscous  forces  at  scales  smaller  than  the  Kolmogoroff  length  (q) 
and  by  inertial  forces  at  scales  larger  than  q where 

n = Cv^/e)^^^  . (1) 


The  mean  rate  of  energy  dissipation  per  unit  volume,  e,  is  given  by 


c E 2 V ejjCy 


where  e^^  is  the  symmetric  rate  of  strain  tensor 


*1J  =2  • 


(2) 


(3) 


The  second  hypothesis  states  that  for  lengths  much  larger  than  q,  the 
statistical  distribution  of  velocl'ty  fluctuations  is  determined  by  c alone. 
This  hypothesis  Implies  that  the  power  spectrum  of  velocity  fluctuations  is 
given  by 
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(J)(k)  = a , L~^  «k  «n~^  (4) 

where  k is  wave  number  by  k = 2TTf/U.  The  range  of  k for  which 

eqn.  (4)  is  valid  is  called  the  "Inertial"  power  law  subrange.  There  is  a 
third  similarity  hypothesis  which  slightly  modifies  the  -5/3  exponent  but 
is  not  of  significance  in  this  paper.  If  one  were  to  make  eqn.  (4)  dimension- 
less using  e and  v one  obtains  the  "universal"  equation  of  all  turbulent  flow 
of 


(()j.(k)  = a (5) 

where  ci  - 0.55  for  a 1-D  spectrum.  The  experimentally  measured  values  of  1-D 
ot  range  from  0.3  to  0.6  with  low  numbers  generally  coming  from  "low"  Reynold's 
number  CRe)  flows  with  the  author  estimating  the  "correct"  value  as  0.55. 

Actual  turbulent  b.l.  have  a finite  Re  and  the  spectral  form  is  illustrated 
in  Figure  1.  The  largest  length  scale  of  the  inertial  subrange  is  the 
boundary  layer  thickness  (6)  and  the  smallest  scale  is  lOq  . At  scales 
smaller  than  lOq,  viscous  forces  become  significant  and  the  spectrum  exponent- 
ially "rolls  off."  An  experimentally  measured  and  normalized  velocity  spectrum 
is  illustrated  in  Figure  2.  A radiated  noise  spectrum  reported  by  Haddle  and 
Skudrzyk  (1969)  is  superimposed  on  the  velocity  spectrum.  The  acoustic  spec- 
trum was  shifted  parallel  to  the  x and  y axis  to  obtain  the  best  fit.  The 
curve  fit  Implies  a frequency  of  46  kHz  as  corresponding  to  n or  0 = c/20TTf  = 50y. 
It  is  the  author's  experience  that  n will  range  from  10  to  lOOy  in  water  so  50y 
is  a reasonable  value. 

The  dilemma  is  illustrated  in  Figure  3 which  illustrates  Haddle  and 
Skudrzyk' s original  three  curves  with  slope  lines  added.  The  radiated  noise 
curves  fit  the  predicted  -5/3  slope  but  the  wall  pressure  spectrum  fits  a 
-11/3  slope.  It  is  suggested  that  the  actual  wall  pressure  spectrum  is 
identical  in  form  to  the  radiated  noise  spectrum  but  the  measured  wall 
spectrum  is  distorted  by  sensor  averaging. 
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Wall  Pressure  Sensor  Correction 

The  data  of  Figure  3 was  obtained  with  D = 2.54  cm  sensors.  At  30  kHz 

L = 2c/TTf  = 3.2  cm  but  L = 2U/TTf  = 0.03  cm.  The  spatial  resolution  of 
a V 

the  sensors  is  adequate  for  the  flow  noise  but  grossly  inadequate  for  the 
wall  pressure.  If  one  knew  what  the  sensors  measured,  one  could  apply  a 
correction  to  the  measured  data  to  get  the  spectral  form  of  the  actual 
phenomenon . 

The  sensor  is  a flat  circular  disk  in  contact  with  a 3-D  eddy.  It  will 
be  assumed  that  the  eddy  has  a circular  cross  section  and  that  the  pressure 
distribution  of  the  eddy  is  given  by  cos  wr  for  -tt/2  1 wr  1 Tr/2.  For  this 
analyses  all  distances  are  normalized  by  the  sensor  radius  and  pressure 
amplitude  by  the  peak  amplitude  of  the  eddy  of  interest.  The  measured 
pressure  is  given  by 

I P dA 

p = — J 

m A TT 


cos  (air)  2'Trr  dr 


(6) 


With  a variable  change  of  u = u>r,  eqn  (6)  becomes 


u cos  udu 


A constraint  on  eqn  (7)  is  that  w £ tt/2  so 


(7) 


If  10  > 7r/2  , integration  limit  = v/2 


(0  < Tr/2  , integration  limit  = (o 


The  integral  of  eqn  (7)  for  (o  > Tr/2  is 

. 2 -1). 

*^m  2 

U) 

and  for  to  < Tr/2 


(8) 


For  (0  * Tr/4,  p = 
m 

falls  6 dB/ octave 


(cos  to-  + to  sin  (0  ~1) 

2 

to 

-0.7  dB  (re  actual  value),  at  to 

_2 

at  higher  frequencies  (f  )• 


(9) 

= Tr/2,  p »r.3.35  dB  and 
m 

The  critical  diameter  or 
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length  of  a sensor  (-3dB)  for  velocity  pressure  measurements  is  • 2U/Trf 

and  for  sound,  ri  = 2C/Trf.  From  Eqn.  (8)  it  is  seen  that  for  sensors  larger 

_2 

than  n > the  measured  signal  is  attenuated  proportional  to  o) 

V • d 


If  one  has  a spectral  distribution  of  f and  measures  it  with  a sensor 
-2  -11/3 

having  sensitivity  of  f , one  will  measure  a spectrum  of  f or  ""-11  dB/octave. 

Referring  back  to  Figure  3,  the  location  labeled  U/nD  at  *■  200  Hz  corresponds  to 

—5/3 

uj  = tt/A  C=  10.7dB)  so  one  could  not  expect  to  get  a f wall  pressure  spectrum, 
the  sensor  was  larger  than  the  largest  eddy. 


Figure  4 illustrates  a wall  pressure  spectrum  measured  by  Foxwell  and 

Grandshaw  (1968)  with  a 0.3175  cm  (1/8  inch)  diameter  sensor.  The  region  from 

-5/3 

100-1500  Hz  approximates  an  f inertial  subrange  and  1.5-2. 5 kHz  approxi- 
mates an  f slope.  The  f and  f lines  are  separated  by  3.3  dB 

at  the  frequency  labeled  2U/itD  which  is  consistent  with  the  above  theory. 

Poxwell  and  Grandshaw  also  made  similar  measurements  with  sensors  2,  4 and  8 
times  this  diameter.  The  curve  for  the  sensor  twice  as  large  had  the  same 
spectral  shape  with  the  break  between  the  two  slopes  at  - 750  Hz  and  at  ~ 350  Hz 
for  x4  sensor. 

Figure  5 illustrates  experimental  data  reported  by  Foxwell  (1966)  taken 
on  a rotating  cylinder.  The  data  was  taken  with  4 sensors,  D = 0.3175  cm  (1/8 
inch),  x2,  x4  and  x8  in  diameter.  The  data  was  corrected  on  the  basis  of  +6  dB/ 
octave  with  an  additional  +6  dB  for  each  doubling  of  the  sensor  size  to  obtain 
the  "corrected"  data.  The  "corrected"  data  of  the  smallest  two  sensors 
collapse  to  a single  line  as  expected  but  the  data  from  the  largest  two  sensors 
fall  short  of  the  -5/3  line.  The  reason  the  largest  sensors  fail  is  that  they 
are  simply  too  large.  After  applying  the  corrections  to  many  experimental  data 
curves  the  writer  has  found  that  the  corrections  are  only  valid  to  40-45  dB. 

The  corrections  are  not  reliable  at  values  in  excess  of  45  dB.  A similar  set 
of  data  by  Foxwell  is  illustrated  in  Figure  6 with  corrections  by  the  author. 

The  "corrected"  data  of  the  larger  two  sensors  falls  further  short  of  the 
-5/3  line  in  Figure  6 than  in  Figure  5.  This  is  consistent  with  the  theory 
since  the  higher  speed  of  Figure  6 produces  smaller  length  scales  which 
exacerbates  the  spatial  averaging  problem.  The  "corrections"  made  in  Figures 
5 and  6 are  relative  corrections  made  to  keep  the  data  on  the  original  graphs. 
The  actual  corrections  would  shift  the  "corrected"  data  vertically  upwards. 

Figure  7 is  a compilation  of  data  of  several  scientists'  wall  pressure 
spectrum  measurements  with  slope  lines  by  the  author.  It's  clearly  evident  that 
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-11/3 

the  data  converges  to  a single  f 
numbers  and  there  is  some  evidence  of 


line  at  the  highest  normalized  wave 
-5/3 

f data  at  the  lowest  wave  numbers. 


Figure  8 Illustrates  data  reported  by  Kirby  (1968)  on  wall  pressures 

-11/3 

spectra  that  Illustrates  that  the  measured  data  Is  proportional  to  f . 

All  of  the  data  points  to  a problem  of  spatial  averaging  on  wall  pressure 

measurements.  Figure  9 Illustrates  some  recent  data  by  Patrick  (1977)  on 

-11/3* 

wall  pressure  spectra  with  spectral  slope  of  f . Figure  10  Illustrates 

some  recent  flow  noise  data  by  Lauchle  (1977)  with  acoustic  signal  proportional 
-5/3 

to  f . The  data  was  obtained  from  a flush  mounted  transducer  In  a region 
of  laminar  boundary  layer.  This  minimizes  the  problem  of  superimposed 
pressures  of  eddies  and  acoustics. 

Conclusions 

1.  The  spectra  of  turbulence,  wall  pressure  and  radiated  flow  noise 
are  similar. 

2.  Measured  spectra  falls  of  6 dB/octave  from  the  actual  spectra 
when  the  data  Is  spatially  averaged. 

3.  A valid  correction  factor  for  large  circular  sensors  Is  given. 


* The  power  spectrum  of  velocity  has  dimensions  of  velocity  squared  per  Hz. 

One  obtains  a pressure  spectrum  by  multiplying  the  velocity  power  spectrum 
by  p/2  to  get  units  of  pressure  per  Hz.  The  power  spectrum  of  pressure  has 
units  of  pressure  squared  per  Hz.  Most  of  the  papers  reviewed  used  pressure 
spectral  density  with  units  of  pressure  per  Hz  which  has  a 1:1  correlation 
to  the  velocity  power  spectrum  density.  The  concept  of  dB  Is  not  very 
appropriate  to  use  In  correlating  velocity  an^  sound.  The  ordinate  labeling 
of  Patrick's  data  has  been  changed  from  dB  (p^/Hz  basis)  to  p/Hz  form  to  be 
consistent  with  the  paper. 
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FIGURE  2.  Spectral  Form  of  Turbulence  Compared  to 
Radiated  Noise 
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FIGURE  3.  Spectra  of  Wall  Pressure  and  Radiated  Noise 
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FIGURE  6.  Spectrum  of  Wall  Pressure  by  Foxwell  at 
U-23.93  m/sec 
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FIGURE  7.  Normalized  Wall  Pressure  Spectra  of  Various 
Scientists 
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FIGURE  9.  Wall  Pressure  Spectral  Shapes  by  Patrick 


FIGURE  10.  Flow  Noise  Measurements  by  Lauchle 
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